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Abstract 

We reconsider the problem of option pricing using historical proba- 
bility distributions. We first discuss how the risk-minimisation scheme 
proposed recently is an adequate starting point under the realistic as- 
sumption that price increments are uncorrelated (but not necessar- 
ily independent) and of arbitrary probability density. We discuss in 
particular how, in the Gaussian limit, the Black-Scholes results are 
recovered, including the fact that the average return of the underly- 
ing stock disappears from the price (and the hedging strategy). We 
compare this theory to real option prices and find these reflect in a sur- 
prisingly accurate way the subtle statistical features of the underlying 
asset fluctuations. 
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1 Introduction 



The famous Black and Scholes option pricing theory has two remarkable fea- 
tures: the hedging strategy eliminates risk entirely, and the option price does 
not depend at all on the average return of the underlying asset |], 0, [3[] . The 
second property means that the option price is not simply the actualized av- 
erage of the future pay-off over the historical probability distribution, which 
obviously would depend on the average return. This is even more striking in 
the case of the Cox-Ross-Rubinstein binomial model j|, ^ where the pricing 
measure is completely unrelated to the actual distribution of returns. This 
has lead to a rather abstract and general framework for derivative pricing, 
where the absence of arbitrage opportunities leads, for models where risk 
can be eliminated completely, to the existence of a 'risk-neutral probability 
measure' (unrelated to the historical one) over which the relevant average 
should be taken to obtain the price of derivatives H H|. It is thus a rather 
common belief that the knowledge of the 'true' probability distribution of 
returns is a useless information to price options. The credence is rather that 
the relevant 'implied' value of the parameters should be obtained from option 
market themselves, and used to price other instruments (for example exotic 
options) @,[7[. 

However, in most models of stock fluctuations, except for very special 
cases (continuous time Brownian motion and binomial, both being very poor 
representation of the reality), risk in option trading cannot be eliminated, 
and strict arbitrage opportunities do not exist, whatever the price of the op- 
tion. That risk cannot be eliminated is furthermore the fundamental reason 
for the very existence of option markets. It would thus be more satisfactory 
to have a theory of options where the true historical behaviour of the under- 
lying asset was used to compute the option price, the hedging strategy, and 
the residual risk. The latter is clearly important to estimate, both for risk 
control purposes, but also because it is reasonable to think that this residual 
risk partly determines the bid-ask spread imposed by market- makers. The 
natural framework for this is the risk minimisation approach developed by 
several authors || ||, [H], [□], [12], [13], [14| , where the optimal trading strategy is 
determined such that the chosen measure of risk (for example the variance of 
the wealth balance) is minimized. The 'theoretical' price is then obtained us- 
ing a fair game argument. Note that in this approach, the option price is not 
unique since it depends on the definition of risk; furthermore, a risk-premium 
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correction to the fair game price can be expected in general. From a theoret- 
ical point of view, this would generally be regarded as a lethal inconsistency. 
From a practical point of view, however, we see this as an advantage: since 
the price ambiguity is a constitutive property of option markets, it is inter- 
esting to understand the origin and size of this ambiguity. In this framework, 
the historical probability distribution determines the 'pricing kernel' to be 
used in the option price formula. We show in detail how, in the Black-Scholes 
limit, the average trend indeed completely disappears from the formula, and 
all the classical results are recovered. For more general models, however, the 
independence of the price on the average return is non trivial. 

The outline of this paper is as follows. We first recall the basic steps 
leading to option prices and optimal hedges for a general process with un- 
corrected (but not necessarily independent) increments, which we present 
in terms of a cumulant expansion to show how the Black- Scholes results are 
obtained in the corresponding Gaussian limit. The first cumulant correction 
provides a theory for the volatility smile in terms of the (maturity dependent) 
kurtosis of the terminal price distribution. We compare this theory to real 
option prices (on a liquid market) and find that these option prices reflect in 
a surprisingly accurate way the subtle statistical features of the underlying 
asset fluctuations Jl5|, [TjJ, in particular the persistent nature of the volatility 
fluctuations. 



2 A risk minimisation theory of option pric- 
ing 

2.1 The global wealth balance 

Let us first write the wealth balance equation corresponding to the writing 
of a European call option. At time t = 0, the writer receives the price of 
the option C[x ,x s ,T], on a certain asset which value is x{t = 0) = x . The 
strike price is x s . Between t = and t = T, the writer trades the underlying 
asset at discrete times t = kr, k = 1,...,N = T/t; his strategy is to hold 
4>k{xk) assets if the price is x(t) = Xk when the time is t = kr. It is easy to 
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show that the change of wealth due to this trading is given by (TT| : 

N-l 

Aiy trading = J2 4kM[xk+i - e^x^-^ (1) 

fc=0 

where r is risk-free rate and t k = kr. At time T = Nt, the writer looses 
the difference — x s if the option is exercized. Thus the complete wealth 
balance reads: 

N-l 

AW = C[x , x s , T]e rT - m a x(x N - x s , 0) + £ (f) k {x k )5x k e r( ?- tk +^ (2) 

where we have introduced the notation: 5xk = [xk+i — e rT Xk] ■ Note that 5xk 
is posterior to the instant k where <pk is determined. Denoting as (...) the 
average over the historical distribution, the average profit is given by: 

N-l 

(AW) = C[x , x s , T}e rT - (max{x N -x a , 0)) + ^ (M^)) (8x k )e r ^ tk ^ (3) 

k=0 

The fair game requirement then fixes C[xo,x s ,T] such that (AW) = 0. 

Although other interesting definitions could be considered [|TTJ , we restrict 
here to the case where the risk is measured as: 

R = (AW 2 ) - (AW) 2 = (AW 2 ) (4) 

The risk R is always greater than or equal to zero and the minimum is ob- 
tained for a certain optimal strategy 0*, determined by a functional deriva- 
tion of (f|) with respect to with respect to <p(x, t). This determines the option 
price through: 

N-l 

C[x ,x s ,T] = e- rT (m a x(x N -x s ,0)) - £ (<i>* k (x k )) (5x k )e~ rt ^ (5) 

Note that since <fi* depends a priori on our choice of the variance as the 
relevant measure of risk, the price of the option is not unique, but reflects 
(among other things) the operator's perception of risk. 
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2.2 The case of zero excess average return 

In this section we consider the case where the average return of the stock over 
the bond, m = (5xk), is zero. This simplifying hypothesis is often justified in 
practice for small maturities, where average return effects are small compared 
to volatilites, and can be treated perturbatively, as shown in section 

Let P(x,T\xo,0)dx be the probability that the asset value is x at time 
T, knowing that it was xq at time 0. When m = 0, Eq. @ then yields an 
option price independent of the trading strategy: 



C[xq, x s , T] m = 0] = e rT (max(xjv— x s , 0)) = e rT dx(x—x s )P(x,T\xo,0) 

Jx s 

(6) 

Note that our assumption that m = means that the average of P(x, T\xq, 0) 
is not at xq, but at the forward price xoe rT . 

In order to proceed with the risk-minimization, we shall assume that the 
price increments Sxk are uncorrelated random variables, such that (5x k 5xe) = 
<? 2 fik,i, where 8 k ^ is the Kronecker symbol. Assuming that a does not depend 
on k is in general not justified, since it amounts to assuming that share 
price follows an additive random process of constant volatility (but with 
an arbitrary distribution for the increments). Actually, real data is often 
closer (for short maturities) to being an additive random process rather than 
a multiplicative one |TT| , an assumption which does introduce an spurious 
positive skew in the price distribution. In reality, however, a depends on k, 
which reflects ARCH-like effects (or time persistent volatility Taking 
this effect into account would lead to more involved calculations, which can 
however still be completed analytically (Tl[. With these approximations in 



mind, the relevant formula for risk is rather simple: 



N-l 



1 " x roo 

(AW 2 ) = (AW 2 ) + a 2 J2 dxP(x,t h \x o ,0)<j>l(x)t 

- 2e r( - T ~ tk +^ J2 dxP(x,t k \x ,0)Mx) 

/ dx'(x' - x s )P{x', T\x, t k ) {5x k ) X)tk ^ x > iT (7) 

where (AW 2 )q is the unhedged = 0) risk associated to the option, and 
{8xk)x,t k ^x',T is the conditional average of Sx k , on the trajectories starting at 
x at time t k and ending at point x' at time T. 
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The optimal trading strategy is obtained by setting JT0], [11], 0: 

d(AW 2 ) 



d(p k (x) 



(8) 



for all k and x. This leads to the following explicit result for the optimal 
hedging strategy: 

e ~r(T-t k+1 ) „oo 

( t>l{ x ) = 2 / dx'(x' - x s )(5x k ) xA ^ x ^ T P(x\T\x 1 t k ) (9) 

(J JX S 

This formula simplifies somewhat when the increments are Gaussian, and 
one finally finds the famous Black-Scholes 'A— hedge': <j)%{x) = dC[x, x s , T — 
tk]/dx. In the non-Gaussian case, however, this simple relation between the 
derivative of the option price and the trading strategy no longer holds (see 
Eq. (0) below). 

Inserting (0) into (0) leads to the following formula for the residual risk: 

N-i r00 

R* = (AW 2 ) -Dry dxP(x, t k \x , O)0f (x)e- r{T -* fc+l) (10) 
k=o Jo 

In general, the left-hand side of (|10D is non-zero; in practice it is even quite 
high - for example, for typical one-month options on liquid markets, V R* 
represents as much as 25% of the option price itself ||11|| . However, in the 
special case where P(x, t\xo, 0) is normal (or log-normal), and in the limit of 
continuous trading, that is, when r — * 0, one can show that the residual risk 
R* actually vanishes, thanks to a somewhat miraculous identity for Gaussian 



integrals ||I0| . Hence, the above formalism matches smoothly with all the 
Black-Scholes results in the limit of a continuous time Brownian (or log- 
Brownian) process, at least when the excess average return of the asset is 
zero. Let us now discuss how these results are changed if the average return 
m = (5x k ) is non zero (but small). 



2.3 Small non-zero average return 

More precisely, we shall consider the case where mN <^ a\^N (N = T/r), 
or, more intuitively, that the average return on the time scale of the option 
is small compared to the typical variations, which is certainly the case for 
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options up to a few monthsQ The global wealth balance then includes the 
term related to the trading strategy, which reads^: 

N-l ,00 

(A^ading) = ™ E / dxP(x, t = kr\x , 0W k (x) (11) 
fc=0 Jo 

The advantage of considering a small average return is that one can do a 
perturbation around the zero average return case, and still use the explicit 
optimal strategy of (^) to lowest order in m. 

Compared to the case m = 0, the option price is changed both because 
P( 0) is biased, and because (AWtrading) must be substracted off from 

Eq. (pp. It is convenient to use the Fourier transform of the probability 
distribution P(z) = /f^ dxP(x, N\xo, 0) exp(zxz) and to expand it in a series 
introducing the cumulants c n . They are defined by 

P(,)=exp[£^], (12) 

n=l U - 

where C\ = mN, c 2 = Na 2 , k = c^/c 2 . is the kurtosis, etc... Applying the 
cumulant expansion to the probability distribution in Eq.(^), we obtain the 
following optimal strategy for m = flTT| : 



1 00 (—lYc (9 n_1 

«(*) = - ^M^^'- T -^ < 13 ' 

Note again that in the Gaussian case (c n = for n > 2), one recovers the 
standard Black and Scholes 'A— hedge'. 

Inserting the optimal strategy flT3| ) into Eq. ([□]) and integrating by parts 
one gets an expansion of the trading term (AWtrading)- Inserted into Eq.(||), 



this gives the following correction to the option price [11 



C[x ,x s ,T;m] = C[x , x s , T; m = 0] (14) 

-iSA*^^ NM '=- + 0(m2) 



1 Typically, m = 5% annual and a — 15% annual. The order of magnitude of the error 
made in neglecting the second order term in m is m 2 N /a 2 ~ 0.1 even for Nr — 1 year. 

2 For the sake of simplicity, we shall set the interest r to zero in the following. See 
for a more complete discussion. 
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where the shorthand Pq stands for the probability distribution where the first 
cumulant m has been set to zero. 

In the Gaussian case, c n = for all n > 3, and one thus sees explicitly 
that C m = Co, at least to first order in m. Actually, one can show that this 
is true to all orders in m in the Gaussian case, which is an alternative way 
to derive the result of Black and Scholes in a Gaussian context jO]] []. 

However, for even distributions with fat tails (c 3 = and c 4 > 0), it is 
easy to see from the above formula that a positive average return m > 
increases the price of out-of-the-money options (x s > xq), and decreases the 
price of in-the-money options (x s < xq). Hence, we see again explicitly that 
the independence of the option price on the average return m, which is one 
of the most important result of Black and Scholes, does not survive for more 
general models of stock fluctuations. 

Note finally that Eq. flT5| ) can also be written as: 

poo 

C[xo, x s , T; m) = / dx(x — x s )Q(x, T\xq, 0) (15) 

with an effective distribution Q defined as: 

m 00 c <9 n_1 
Q(x, T\x , 0) = P (x, T\x , 0) - - E 7 ^Tua n I P o( x > N \ x o> °) ( 16 ) 

c 2 n=3 \P>— J-J ! ox 

The integral over x of Q is one, but Q is not a priori positive everywhere. 
This means that for a certain family of pay- offs, the fair price of the option 
flTBT) may be negative in the absence of risk-premium. From a practical point 
of view, however, this requires rather absurd values for the average return 
and for the strike price, which in turn would lead to a large residual risk. 

The pseudo-distribution Q generalizes the 'risk neutral probability' usu- 
ally discussed in the context of the Black-Scholes theory, and also has the 
property that the excess average return (the integral of (x — x )Q over x) is 
zero, as can easily be seen by inspection from flT6|) . In fact, one can derive 
a general formula for Q without any restriction on m of r [fTI^ , 13 1, and the 



effective distribution still has the properties that it is normalized to one and 
has zero average excess return. 



3 This result is obvious on the framework of Ito's calculus, which is only valid for all 
Gaussian processes (including the log-normal) in the continuum time limit. 
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3 Volatility smile and implied kurtosis 



In the case where the market fluctuations are moderately non-Gaussian, as 
is the case for liquid markets, one might expect that the first terms in the cu- 
mulant expansion around the Black-Scholes formula are sufficient to account 
for real option prices. If one only retains the leading order correction which 
is (for symmetric fluctuations) proportional to the kurtosis k, one finds that 
the price of options C(xq,x s ,T) can be written as a Gaussian Black-Scholes 
formulafj, but with a modified value of the volatility a, which becomes price 
and maturity dependent |L7] : 



v imp (x s ,T) = a 



24 I a 2 T 



(17) 



The volatility o"i mp is called the implied volatility by the market operators, 
who use the standard Black-Scholes formula to price options, but with a 
value of the volatility which they estimate intuitively, and which turns out 
to depend on the exercice price in a roughly parabolic manner, as indeed 
suggested by Eq. (|rT|). This is the famous 'volatility smile'. Eq. (|I~7l) 
furthermore shows that the curvature of the smile is directly related to the 
kurtosis k t of the underlying statistical process on the scale of the maturity 
T = Nr. We have tested this prediction by directly comparing the 'implied 



kurtosis' ||18|| , obtained by extracting from real option prices the volatility a 
and the curvature of the implied volatility smile, to the historical value of the 
volatility and of the kurtosis k^. We have mostly studied short maturities 
(up to two months) options on futures, for which the interest rate can be set 
to zero. We also restrict to liquid markets (such as the bund option market) 
where (i) non Gaussian effects are not too strong, and (ii) risk-premiums are 
expected to be small, and thus where a comparison with the fair price is 
meaningful. 

We have found the following results. The implied volatility turns out to be 
highly correlated with a short time filter of the historical volatility: see Fig. 1. 
Fig. 2 shows the comparison between the implied and historical kurtosis, with 



4 Note that the operators rather use the more standard log-normal Black-Scholes for- 
mula, which, as noted above, induces a spurious positive skew not present in real data 
(at least for short maturities). In order to correct for this skew, the log- normal volatility 
smile is then negatively skewed. A more symmetric smile is observed if one talks in terms 
of a Gaussian volatility, which is what we adopt in the following. 
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Figure 1: Comparison between the historical volatility of the bund (measured 
from high frequency data and filtered over the past five days), and the implied 
volatility, extracted from the option prices through formula (|l7l). 
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no further ajustable parameters. Note that the historical kurtosis decays more 
slowly than iV _1 , which would be expected for a process with independent, 
identically distributed increments. This anomalously slow decay is directly 
related to volatility persistence effects [17], [Tl| . 

It is interesting to note that the kurtosis correction to the optimal strategy 
does not coincide with the market practice of using the implied volatility in 
the Black-Scholes A— hedge. However, since the risk is minimum for <fi = </>*, 
this means that the increase of risk due to a small error 5<p in the strategy is 
only of order 5<p 2 , and thus often quite small in practical applications. 

The remarkable agreement between the implied and historical value of 
the parameters (which we have also found on a variety of other assets), and 
the fact that they evolve similarly with maturity, shows that the market 
as a whole is able to correct (by trial and errors) the inadequacies of the 
Black-Scholes formula, and to encode in a satisfactory way both the fact 
that the distribution has a positive kurtosis, and that this kurtosis decays 
with maturity in an anomalous fashion due to volatility persistence effects. 



4 Conclusion 

In our opinion, mathematical finance in the past decades has overfocused on 
the concept of arbitrage free pricing, which relies on very specific models (or 
instruments) where risk can be eliminated completely. This leads to a re- 
markably elegant and consistent formalism, where derivative pricing amounts 
to determining the risk-neutral probability measure, which in general does 
not coincide with the historical measure. In doing so, however, many impor- 
tant and subtle features are swept under the rug, in particular the amplitude 
of the residual risk. Furthermore, the fact that the risk-neutral and histor- 
ical probabilities need not be the same is often an excuse for not worrying 
when the parameters of a specific model deduced from derivative markets are 
very different from historical ones. This is particularly obvious in the case of 
interest rates |I9|. In our mind, this rather reflects that an important effect 



has been left out of the models, which in the case of interest rates is a risk 



premium effect |19| . We believe that a more versatile (although less elegant 
from a mathematical point of view) theory of derivative pricing, such as the 
one discussed above, allows one to use in a consistent and fruitful way the 
empirical data on the underlying asset to price, hedge, and control the risk 



11 



10 



V 1 - 




i Historical Kp (future) 

) Implied K. mp (option) 
Exponential vol. correlations 
Power-law vol. correlations 



10 



100 
N 



1000 



Figure 2: Plot (in log-log coordinates) of the average implied kurtosis Ki mp (de- 
termined by fitting the implied volatility for a fixed maturity by a parabola) and 
of the empirical kurtosis kn (determined directly from the historical movements 
of the BUND contract), as a function of the reduced time scale iV = T/r, r = 30 
minutes. All transactions of options on the BUND future from 1993 to 1995 were 
analyzed along with 5 minute tick data of the bund future for the same period. We 
show for comparison a fit with kn — N~ - 6 (dark line). A fit with an exponentially 
decaying volatility correlation function is however also acceptable (dotted line). 
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the corresponding derivative security. Extension of these ideas to interest 
rate derivatives is underway. 
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